Age-dependent transient shear banding in soft glasses 
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We study numerically the formation of long-lived transient shear bands during shear startup within 
two models of soft glasses (a simple fluidity model and an adapted 'soft glassy rheology' model). The 
degree and duration of banding depends strongly on the applied shear rate, and on sample age before 
shearing. In both models the ultimate steady flow state is homogeneous at all shear rates, consistent 
with the underlying constitutive curve being monotonic. However, particularly in the SGR case, 
the transient bands can be extremely long lived. The banding instability is neither 'purely viscous' 
nor 'purely elastic' in origin, but is closely associated with stress overshoot in startup flow. 
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Many soft materials show shear banding: the separa- 
tion into layers of different shear rate under imposed flow. 
Examples include wormlike micelle solutions [![; granu- 
lar matter Q; star polymers H; and 'soft glasses' such 
as gels, pastes and emulsions (4J. Because shear banding 
dramatically alters the stress response at fixed shear rate 
or vice versa, it is central to the rheological control of 
such materials, whose applications range from foodstuffs 
and pharmaceuticals to paints and well-bore fluids. For 
materials that are nonergodic at rest, unexpected com- 
plexity can arise when one band is not flowing and thus 
subject to aging [3H6fl, while the other is continuously 
rejuvenated by flow [7j. This interplay between glassy 
dynamics and shear banding brings together two major 
current areas of nonequilibrium physics research. 

A steady mean shear rate is imposed by choosing time- 
independent wall velocities in a rheometric device. In 
most reported cases, the resulting shear bands then per- 
sist indefinitely In many instances, shear 
banding is attributable to a nonmonotonic steady-state 
constitutive curve £(7), where £ is the shear stress and 
7 the shear rate in the homogeneous material; regions 
with cE/dj < are mechanically unstable [l|. This is 
analogous to, but distinct from, a similar instability in 
nonlinear elastic solids, where mechanical instability like- 
wise arises if the stress is a decreasing function of strain 
{dH/d/y < 0); see 13|. For simplicity we refer to these as 
'viscous' and 'elastic' banding scenarios respectively. 

In this Letter we use the case of soft glasses to explore 
theoretically a distinctive, third scenario. This is where 
shear banding arises on startup of steady shearing, before 
eventually relaxing to an unbanded steady state. The re- 
sulting transient bands can be extremely long-lived and 
so might well be mistaken for true steady-state ones (see, 
e.g. [14{). We show that such bands can arise in systems 
showing neither a viscous nor an elastic instability, such 
as a model combining essentially linear elasticity with 
a near-trivial constitutive curve, £(7) = A + Bj. Our 
results may thus prove highly relevant (alongside other 



factors [15|, |16|) to a number of cases where apparently 



steady shear banding is seen, even though the consti- 
tutive curve is pre dicted , by well founded theories, to 

They are also relevant to re- 



remam monotonic 
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cent experiments where long-lived transient shear bands 



are directly reported [21 



We argue that transient banding behavior can arise 
generically in systems where the stress response £(£,7) 
to startup of steady shear shows, as a function of time, 
a significant overshoot before falling to the steady-state 
limiting value £(00,7) = £(7). Indeed, the simplest case 
is where the overshoot is created purely by nonlinear elas- 
ticity without relaxation or plastic flow. In this case, 
7(f) = jt is an elastic strain, so that any region where 
9t£(£, 7) < implies the onset of the elastic instabil- 
ity referred to already. In reality however, soft glasses 
(in contrast to, e.g., elastomers 13[) have a very limited 
elastic deformation regime before they yield irreversibly; 
the stress maximum in such cases is the result of a con- 
test between the growth of elastic stress and its decay by 
plastic rearrangement towards the eventual steady-state 
limit. Accordingly 7(4) is not an elastic strain within 
the region where <9j£ < 0, and there can be no direct 
mapping onto an elastic banding scenario. 

Below we show that transient shear banding does 
nonetheless arise, whenever the stress overshoot becomes 
large, both in a fluidity model of soft glasses (chosen 
for its tractability) , and in a mesoscopic model, adapted 
from that of [22[ which is known to capture well some sub- 
tler physics of these systems. In soft glasses the overshoot 
can be varied in height (without also varying £(7)) purely 
by changing the age of the system [22| , making such ma- 
terials ideal testing grounds for experiments and theory 
on transient banding. Emerging as it does from two in- 
dependent models, transient shear banding should be a 
generic feature in the rheology of well-aged soft glasses - 
a fact not appreciated previously. Moreover, stress over- 
shoots often arise in other types of soft matter, and in 
many of these, similar mechanisms may be at work. 



Our first model is an empirical fluid- 
lj| 2(j, in which a single 



Fluidity model: 
ity model, along the lines of 
structural relaxation time r regresses continuously to- 
wards a steady-state value determined by the local shear 
rate. A finite diffusivity for r is also introduced, so as to 
prevent band interfaces from becoming infinitely sharp. 
This model shows transient banding in the overshoot re- 
gion, with a strong dependence on system age, which sets 
the initial value of r and thus the overshoot height. How- 
ever the lifetime of the bands apparently remains mod- 
est, regardless of the system age before shearing begins, 
in contrast to the adapted SGR model considered later. 

To set up the fluidity model we decompose the total 
shear stress £ into a viscoelastic stress a arising from the 
glassy degrees of freedom and a purely viscous part: S = 
(7+777. We assume translational invariance in the x (flow) 
and z (vorticity) directions but allow nonuniformity to 
develop in the flow gradient direction, y. Force balance at 
zero Reynolds number (neglecting inertia) then requires 
£ to be independent of y. We suppose a Maxwell-type 
constitutive equation for the viscoelastic stress 



dt<J = Gj — a jr 



(1) 



where G is an elastic modulus and r is the structural 
relaxation time (inverse fluidity) with its own dynamics: 



d t T = l~T/f(j)+l 2 d 2 y l 



(2) 



In the absence of flow this represents simple aging (t = t) 
but with flow present, aging is cut off at the inverse strain 
rate. (This is characteristic of strain-induced plasticity 
0.) Here l a is a mesoscopic length (which could depend 
on 7 without affecting our presentation) describing the 
tendency for the relaxation time of a mesoscopic region 
to equalise with those of its neighbors. 

We choose the steady-state relaxation time as f = To + 
A/|7j, so that in steady state £ = (7 + 777 with a = Gjf = 
GX + Gto7, has a yield stress beyond which it is trivially 
monotonic (a Bingham fluid). We consider flow between 
infinite flat parallel plates at y = 0, L y , and rescale strain, 
stress, time and length so that A = G = tq = L y = 1. 
The steady-state constitutive curve is then 



£(7) = ! + (!+ 77)7 



(3) 



where for simplicity we now set 77 = 0.05. 

We study a shear startup protocol defined as follows. 
First imagine preparing the sample by a deep quench at 
time t — 0, which we assume results in a fully rejuvenated 
initial state with r(y,t = 0) = l,a(y,t = 0) = across 
the whole sample. Next we allow the sample to age at 
rest (so t = 1) until a time i w , before setting the upper 
plate moving along x with constant speed jL y . This 
defines the average imposed shear rate 7 = J j(y, t)dy. 

Transient stress curves £(£,7) for shear startup in a 
homogeneous system are shown for several different t w 
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FIG. 1: (a) Homogeneous startup transient in fluidity model 
at applied shear rate 7 = 0.1 for waiting times t w — 
10°, 10 2 , 10 4 , 10 6 , 10 8 (t w increasing with increasing size of 
overshoot) showing stable region (solid), unstable (dashed). 
Dotted curve: actual startup curve arising from inhomoge- 
neous flow with lo = 0.01 and S = 0.01. (b) Counterpart 
in modified SGR model for x — 0.3,7 = 0.1, iu = 0.05, n = 
50,m = 300, and waiting times t w = 10°, 10 2 , 10 4 , 10 6 . Solid 
lines: banding disallowed. Dotted: banding allowed. Inset: 
approach to steady state at later times in the banding case. 



in Fig. QJi. There is an overshoot that depends strongly 
on the age t w of the sample; one may show that it occurs 
at a strain 7 = j(t — t w ) given by 7 Q exp(7 G ) = 7/j w . To 
within a logarithmic correction this gives ~f = log(7f w ), 
which we use for convenience below in preference to the 
full implicit expression. The response of the sample prior 
to (but not beyond) the stress maximum is almost elastic, 
so the peak stress obeys £ ~ 7 Q . 

To gain intuition for the consequences of the over- 
shoot, consider now a thought experiment in which one 
tracks this homogeneous stress transient for several iden- 
tical sample replicas, each subject to a different shear 
rate. At any fixed time interval t — t w this defines an 
'instantaneous constitutive curve' £(7), with the tilde 
denoting a dependence on t — t w and t w , suppressed for 
clarity of notation. As shown in Fig. [2ja), each such 
curve is non-monotonic over a typical shear rate window 
7 < Jo/(t — t; w ), with monotonicity being restored at 
higher shear rates. This can be understood as follows. 
At any fixed time interval t — t w , those replicas at high 
shear rate j(t — t w ) > 70 are expected to have reached 
steady state, with stresses obeying the constitutive curve, 
Eq.©. In contrast those for which j(t — t w ) < 70 will 
still be on the elastic branch of the stress transient, with 
a stress a = j(t — t w ) 3> £(7). With this transient non- 
monotonicity in mind, we performed an instantaneous 
linear stability analysis about the evolving homogeneous 
flow, determining for each t — t w an instantaneous eigen- 
value whose positivity indicates instability to the onset 
of banding. The unstable windows are shown by dashed 
lines in Fig. QJ,, and correspond broadly to the regions 
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FIG. 2: (a) For chosen t w = 10 6 , plots of the instanta- 
neous homogeneous constitutive curves during start up for 
t-t w = 40, 80, 120, 160, 200 with regions of instability dashed, 
(b) Strength of banding (twx — 7mm) against t — t w for 
7" = 0.1, and tw = 10 4 , 10 6 , 10 s . (c) For the same 7" = 0.1, 
and t w = 10 4 (solid), 10 6 (dotted), 10 8 (dashed), plots in 
the t — t w , 7 plane showing the upper (thick lines) and lower 
(thin lines) boundary curves of the window in which band- 
ing is significant (threshold set by ("y m ax — 7mm)/7 > 0.01). 
(d) Snapshots of banded profiles for 7 = 0.1, t w = 10 6 , and 
t — t w = 5, 120, 130, 135, 225 showing onset, development and 
eventual collapse of the bands. / = 0.01, 8 — 0.01 in (b-d). 



of negative slope in £(7). However, they do not do so 
exactly; this cannot be viewed solely as a 'viscous' insta- 
bility of the instantaneous constitutive curve. Nor does 
the window exactly correspond to that of negative slope 
in E(i, 7), as it would for an 'elastic' instability. Al- 
though clearly associated with a large stress overshoot, 
the banding scenario found here is thus distinct from ei- 
ther the viscous or the elastic one - a view reinforced by 
the linearity of the step strain response in ([1]), and the 
trivial monotonicity in ([3]), for the underlying model. 

We now turn to study the full heterogeneous dynamics 
of the model in this shear startup protocol, allowing spa- 
tial variations in the flow gradient direction y. In each 
run we add a small perturbation a(y, t = 0) = 6 cos(7ry) 
where 5 <C 1, in order to trigger any banding. In each 
run we tracked, as a function of shearing time t — t w , 
the degree of shear banding in the sample, as measured 
at any instant by the difference 7 ma x — 7min between the 
maximum and minimum shear rates present in the cell, as 
shown in Fig. [2jb) for several different waiting times t w . 
Regions of high 7 max — 7 m i n broadly match up with the 
regions of negative slope in Fig.[5Ja), consistent with the 
idea that this instability in the instantaneous constitu- 
tive curve indeed triggers banding. Likewise in Fig. [2jc) 
we delineate the temporal windows, for various i w and as 
a function of shear rate 7, within which banding is sig- 
nificant, as defined by a criterion (7 max — 7min)/7 > 0.01. 
The shear stress that arises in these transiently banded 



flows is of course different from the predictions of the 
homogeneous model; both are shown in Fig. as a 
function of t — i w , at a fixed shear rate 7 = 0.1 for vari- 
ous waiting times. Snapshots of the transient bands for 
one particular startup are shown in Fig. [5] (d). Note 
that the minimum shear rate, which arises within the low 
shear-rate band, can be negative during transient band- 
ing. This is because the low shear-rate band is well-aged, 
has been subjected only to a weak flow, and therefore 
remains essentially an elastic slab. However the stress 
X on this slab decreases post-overshoot: like any elastic 
solid being unloaded, it shears backwards. 

Adapted SGR model: Our second model is a spatially 
resolved adaptation of the SGR model @, in which 
a spectrum of jump rates describes hopping over strain- 
modulated local rearrangement barriers at an effective 
noise temperature x = T/T g ; elastic strain builds up 
and then is released in these plastic jump events. Fol- 
lowing Ref. [6], where full details can be found, nu- 
merically we take j — l...m SGR elements on each of 
i = l...n streamlines, corresponding to y = 0...1, with 
periodic boundary conditions. The stress on streamline 
i is Oi = {k/m)^2ijlij, with elastic constant k = 1. A 
waiting time Monte Carlo algorithm is used to choose 
stochastically the next element to jump. Supposing the 
jump occurs at element ij when its local strain is i = I, 
force balance is then imposed by updating all elements on 
the same streamline as I — > 1 + l/m, and further updat- 
ing all elements throughout the system as I — s> I — l/mn. 
In contrast to Ref. 6|, here we take the noise tempera- 
ture x as constant, but include instead a stochastic jump- 
induced straining of elements on neighboring streamlines 
after each jump event: further adjusting the strain of 
three randomly chosen elements on each adjacent stream- 
line i± 1 by lw{— 1, +2, —1). This creates a diffusive cou- 
pling of the dynamics on different streamlines (different 
y), analogous to that in Eq.([5]), with a strength set by w. 
It also mildly alters the constitutive curve from that of 
pure SGR, without losing monotonicity. 

The evolution of stress with time in startup flow is 
shown in Fig. [TJj, both for the case where homogeneous 
flow is imposed (solid lines), and for the full dynamics 
allowing banding (dotted lines). The trends are quite 
similar to those from the scalar model: for older sam- 
ples, the dotted line departs from the solid one just after 
the overshoot. The initial effect of transient banding, as 
before, is always to decrease the stress to values below 
that of the homogeneous system at the same time point. 

There is however a new feature in the SGR model, 
not seen in the fluidity model. For the oldest samples 
(t w ~ 10 6 , measured in units of the microscopic attempt 
rate for jumps) the time scale for the stress signal to de- 
cay to the limiting value (corresponding to homogeneous 
flow) is inordinately long: indeed, this decay can require 
strains jt of order thousands, as opposed to the order- 
unity values seen in the fluidity model. Since strain re- 
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FIG. 3: Shear rate profiles corresponding to slowly decaying 
stress signal at t w = 10 6 in Fig. [lb averaged over time win- 
dows 250 - 500; 500 - 1000; 1000 - 1500; 1500 - 2000; 2000 - 
2500; 2500 - 3000 (top to bottom at y = 0.5). 



in other classes of viscoelastic soft matter at high enough 
shear rates, including non-aging systems such as entan- 
gled polymers. Indeed, some studies of polymeric mate- 
rials and models have started to explore the connection 
between this and transient banding (l7l[l8j|. Nonetheless, 
we hope that further experimental and theoretical work 
on soft glasses will help elucidate this connection in more 
general terms. Arguably these represent ideal materials 
with which to explore the problem, because the size of 
the overshoot can be varied without any change to the 
final steady state of the system. In all systems both the 
final state and the overshoot depend on the imposed flow 
rate and on sample composition, temperature etc.; but in 
glasses one can also vary the system age t w which affects 
only the overshoot, and not the steady state. 



juvenates the material, this is possible only because, for 
very old samples, the strain rate in the low-shear band 
remains extremely small compared to 7. This behavior, 
clearly visible in the banding profiles presented in Fig. |3l 
is consistent with having an age-dependent static yield 
stress that can lie well above the homogeneous steady 
state stress £(7) . In principle the banded state might 
then persist indefinitely, but in our model the low shear 
band seemingly is eroded slowly by the spreading of the 
fast band, perhaps as a result of the diffusive nonlocality 
in the jump dynamics. 

Discussion: Aside from its longer-lived transients, the 
adapted SGR model shows broadly similar phenomenol- 
ogy to the fluidity model, supporting our contention 
that the transient shear-banding scenario reported here 
is generic for soft glasses. As previously stated, the soft 
glass case is one where the stress overshoot stems from 
a competition between (essentially linear) elasticity and 
plastic relaxation: there is no facile mapping onto ei- 
ther an elastic instability (dS/dj < 0) or a viscous one 
(dYi/drf < 0). Despite the intermediate character of the 
instability towards transient bands, we showed for the flu- 
idity model that their onset is closely correlated with the 
occurrence a negative slope on the instantaneous stress 
versus strain rate curve (dYl/d^t < 0). However, lin- 
earization about the time-dependent solution for a ho- 
mogeneous flow showed this to be a qualitative rather 
than an exact correspondence. 

The scenario of transient bands we have developed 
for soft glasses may interact in a complicated manner 
with various mechanisms previously presented to explain 
ste ady -state shear bands in the same class of materials 
@, [T(|. However, in some cases bands of long but fi- 
nite duration were already observed in startup flows (2l| 
and the physics we have explored in this letter may be 
enough, on its own, to explain some of these. Even if not, 
the presence of a generic connection between transient 
banding and stress overshoots should not be overlooked 
in future work on these materials. 

Stress overshoots in startup flows are often also seen 
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